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Abstract 

A local perturbation theory for the spectral analysis of the Schrodinger operator 
with two periodic potentials whose periods are commensurable has been constructed. 
It has been shown that the perturbation of the periodic ID Hamiltonian by an ad- 
ditional small periodic potential leads to the following spectral deformation: all gaps 
in the spectrum of the unperturbed periodic Hamiltonian bear shifts while any band 
splits by arising additional gaps into a set of smaller spectral bands. The spectral 
shift, the position of additional gaps and their widths have been calculated explic- 
itly. The applications to the operational regime of a nanoelectronic device based on 
Mott-Peierls stimulated transition have also been discussed. 

1 Introduction 

In series of papers ^ [H El HJ 03 El a new type of nanoelectronic device based on the phe- 
nomenon of the stimulated Mott-Peierls transition (SMPT) has been studied in one elec- 
tron and one-mode approximation. The spectral properties and the operational regimes 
of the device have been calculated for the ballistic electrons moving along the quantum 
wire (QW) treated as the conducting electrode of the device. The effect of QW-structure 
has been taken into account by the effective mass m* in the Schrodinger equation 

^ d 2 ^ 

(1) 



2m* dx 2 



under assumption that normally the conducting electrode is in metallic state. When 
a nontrivial operating voltage has been applied to the governing electrodes distributed 
periodically along the QW it has been shown that the effect of the governing electrodes 
can been reduced to an effective periodic potential V e ff in the Schrodinger operator (JJ): 

h 2 d 2 ^ 

The spectrum of the operator © consists of spectral bands separated by gaps. The device 
is in conducting state if the Fermi level of the conducting electrode gets into one of the 
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artificially created spectral bands. On the contrary, if the Fermi level gets into a gap 
between the spectral bands the conductivity vanishes and thus the conducting electrode is 
effectively transformed into the dielectric state. Summarizing what has been said one can 
conclude that the formation of the energy gaps in the continuous spectrum of the operator 
(0 and the control of the parameters of V e f f may be considered as the physical base of a 
new type of nanoelectronic device. 

For a more detailed analysis of spectral properties and physical characteristics (say, 
conductance) of the device described above deviations from the ideal behavior (one- 
electron, one-mode ballistic regime at zero temperature) should be considered because 
the effect of electron-phonon interaction, electron-electron interaction, stimulated super- 
lattice structure of governing electrode , etc., may become important for real device. In 
particular, if the electron concentration is large enough, electrons may be effected by a 
atomic lattice structure of a QW that could change the behavior of the conductance. 
Also if an artificial lattice (periodic) structure could be fabricated for quantum wire then 
two periodic structures (atomic and artificial) could change the spectral properties of the 
device as well as its operational regimes. 

In the present paper we study the spectral properties of the controlled QW which is 
described by ID quantum Hamiltonian 

d 2 

H = + Qo{x) + Sgx(x) (3) 

with the periodic potentials qo(x) and q\(x) playing the role of the atomic potential of 
QW in one-electron approximation (q ) and the effective potential (qi) generated by the 
periodic system of governing electrodes. Everywhere below we set h/2m = 1 where m 
is the electron mass and in these units the parameter 5 is proportional to the ratio of 
go and q\ amplitudes. We assume that the periods a and b of qo and q\ respectively are 
commensurable and we shall call such QW as ID super lattice. We consider 5 as a small 
parameter which is in agreement with strength of the atomic potential in comparison with 
the amplitude of the V e ff calculated in To study the spectrum of the Hamiltonian 
and its eigenfunctions we develop a special local perturbation theory which allows to 
calculate all spectral ingredients of H not only for regular points but also for singular ones. 
The latter as it is shown in the paper coincide with the ends of spectrum of the operator 
H : 

d 2 

H = + Qo(x) (4) 

and with some set of points inside the bands of Hq. This set is described in the frame 
of our approach explicitly. We prove that the main spectral effect of the perturbation is 
the following deformation of <j{H$): all gaps of H$ bear shifts while any spectral band 
of Hq splits by the arising additional gaps into a number of smaller spectral bands of 
H. It is shown that spectral shifts, the locations of additional gaps and their widths can 
be calculated in arbitrary order of the perturbation theory explicitly. In this respect the 
spectral analysis presented here can be considered as the sound mathematical background 
for numerical calculations of real nanoelectronic devices fabricated on the basis of narrow- 
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gap semiconductors 0. The calculations of the spectral characteristics and operational 
regimes for real SMPT-device will be done elsewhere. 

2 The Perturbation Theory at Regular Points 

We consider the perturbation of the Schrodinger operator with a periodic potential by 
another periodic potential with a small amplitude. So the Hamiltonian has the form 

H = -dl + q (x) + Sqi{x), x G R, 

where qo(x + a) = qo(x) and q\(x + b) = q\{x) are two periodic functions, the parameter 
5 is considered to be small. We assume that the periods a and b are commensurable, i.e. 
b/a = m/n where m and n are integers. Hence, the perturbed Hamiltonian H remains 
to be the operator with a periodic potential and c = bn = am is the common period of 
functions qo and q\. 

In what follows we need some well known facts |5] from the theory of the periodic 
Schrodinger operators H per . The spectrum of such operators consists of spectral bands 
separated by gaps. In the general case there are infinite number of bands and gaps. If 
the period of the potential is A then the so-called Bloch functions of H per are solutions of 
equation H per ^f = E^> which obey the quasiperiodic boundary conditions on the period 
i &±(x + A,k) = ^±(x,k). The quasimomentum k is linked with the energy E by 

the dispersion function E = E(k). Under this mapping to the N — th band of H per there 
corresponds the Brilloin zone [ir(N — 1)/A, ttN/A] of the values of quasimomentum. To the 
N-th gap there corresponds the cut [irN/A, ttN/A + i^n] drawn on the complex plane of 
quasimomentum. The so-called Lyapounov function F{E) = cos kA is the entire function 
of E, \F(E)\ < 1 in bands and \F(E)\ > 1 in gaps. At the end of bands ^E(k) = 
and at the gap points E^ corresponding to the tops ttN/A + ijN of cuts the relation 
^F(E) = is valid. 

We assume that the location of bands of the unperturbed operator Hq = —d%. + 
qo(x), its Bloch functions and the dispersion relation are known. Our aim is to study the 
corresponding quantities for the perturbed operator H and to describe its spectrum in 
comparison with the spectrum of the unperturbed operator. 

Let us present the Bloch solutions ^(x,p) of the perturbed spectral problem 

= EV (5) 

in the form 

*(x,p) = e^ p -^ x ^(x,p), (6) 

where E = E(p) and p = p(E) are the dispersion function and the quasimomentum 
corresponding to the operator H, A = \{k) and k = k(X) are the analogous objects for 
the unperturbed problem 

tfo*o = A* , H = -d 2 x + q (x). (7) 
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Due to Floquet theorem the function e lkx &(x,p) is the periodic one with a period c. 
Hence the function <&(x,p) is quasiperiodic and 

$(x + c,p) = e ikc $(x,p). (8) 

To fix the Bloch solution ^(x,p) and hence the function $>{x,p) unambiguously we 
normalize it as follows 

*(0,p) = 1, *(c,p) = e ipc . (9) 

Obviously, 

$(0,p) = 1, $(c,p) = e ifcc . 
On use of eqs.© and © one can see that the function $(x,p) satisfies the equation 

-d 2 x - 2i(p - k)d x + (p- kf + q (x) + 5q x {x) - e] $(x,p) = 0. (10) 

This equation is to be solved in the class of functions with the boundary conditions (jSJ). 
Let us seek the solution of eq. (fTU|) in the form 



= * (x,k) + Y,8 n * n (x,k) (11) 



n>l 

assuming that the perturbed quasimomentum p and the dispersion function E = E(p) are 
also expanded in series 

p = k + $ n »n(k) (12) 
n>0 

and 

E{p) = \{k) + S n ^n{k), (13) 

n>0 

where k and X(k) are the quasimomentum and the dispersion function corresponding to 
the unperturbed problem (JJJ). 

Since the unperturbed Bloch function ^oix, k) satisfies the condition ^q{x + a,k) = 
e ika ^>o(x,k) and hence ^>o( x + c,k) = e tkc ^> o(x , k) , in order to provide the boundary 
condition © one has to seek the correction terms under the same conditions, i.e. 

^ n (x + c,k) =e ikc $> n (x,k), n> 1. (14) 

The normalization conditions 

* (0,fc) = l, V (a,k) = e ika , (15) 

*„(0,fc) = ^ n (c,k) =0, n> 1 (16) 

provide the property © of the exact Bloch solution. 

In what follows we propose a procedure of calculation of all correction terms ^ n (x, k), 
fj, n (k), \ n (k) and define the values of unperturbed quasimomentum k at which this proce- 
dure fails. The corresponding energy points X s = X(k s ) we shall call as singular points of 
perturbation series constructed here. The complement of a set of singular points on real 
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axis we shall call as a set of regular points. In the next Section we construct the local 
perturbation series which are valid at the above singular points. 

The procedure of calculating the correction terms ^ n (x,k), n n (k) and \ n (k) includes 
the following steps: 

1. For each n > 1 we find the linear link between the terms \ n (k) and /i n (k) in 
expansions (|12j) and To this end we use the fact that these expansions are not 
independent and are linked with each other by means of equation 

T(E) = cospc (17) 

where F(E) is the Lyapounov function of the perturbed operator H. 

2. For each n > 1 we find the second linear link between X n (k) and /J, n (k) which 
is necessary to define X n (k) and fi n (k) unambiguously. This second link is obtained by 
examining the solvability condition of the inhomogeneous equation for $> n (x,k) in the 
class of functions with boundary conditions (|14j>. 

3. Finally we find the solution fy n (x,k) which satisfies the boundary conditions Q14j) 
and (|Tnj). 

All technical details of the realization of the above three steps are contained in the 
proofs of statements formulated below. 

Proposition 1. The Lyapounov function J~(E) associated with the operator H in 
the leading order of 5 can be expressed in terms of the unperturbed quasimomentum k as 
follows ^(E) = cos kc + 0(5). 

Proof. Let us insert the expansion IJ13JI into the r.h.s. of the dispersion relation (|17|) 
and expand the r.h.s. in the Taylor series at the vicinity of k. If we restrict the Taylor 
series to the leading order we obtain the statement of proposition. Q.E.D. 

Remark . It should be noted that the leading term cos kc of the Lyapounov function 
^(E) does not coincide with the Lyapounov function i<b(A) = cos ka of the unperturbed 
operator Hq. This fact has crucial consequences for the spectrum of the perturbed operator 
which is to be discussed below. 

To formulate the next statement we need the following notations: i ? (A)|^ = ^( / t) = cos kc, 
F> = AF(A), Fin) = J^ (A); k (k) = £ A(fc). 

Proposition 2. The correction terms of series and mfy) for E(v) andp are linked 
by the recursion relations 

X(k)n n (k) - \ n (k) = T n {\i,fii,i < n - 1} (18) 

where T\ = const and T n {•,•}, n > 1, are some functions of Aj, < n — 1, which can be 
calculated explicitly for arbitrary n. In particular, the first T n , n = 1,2, 3,have the form 

Ti = 0, 

r a (Ai,m) = i(^A?-p^ 

\(k)c 2 nf F" . x F"' F 2 
T 2 (\i,m, A 2 ,/i 2 ) = + -pf AiA 2 + — X 1 + — 
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Proof. To obtain the above relations one has to insert expansions (|13j) and ()12|) for E 
and p into the dispersion relation (|17jl. to expand the l.h.s into the Taylor series at the 
vicinity of A, to expand the r.h.s. into the Taylor series at the vicinity of k and to equate 
coefficients at each power of 5. Finally one has to take into account that in the leading 
order F{E) = F(X). Q.E.D. 

The ra-th relation in (|18|) gives one linear link between two quantities X n {k) and fJ> n (k). 
To define them unambiguously one needs the second link between them. To obtain it 
let us consider the solvability condition of the equation defining the n-th correction term 
fy n (x,k) of the Bloch function. To write down this equation let us insert the expansions 
(|12j) and (|13|) into the exact equation (jlUj) and equate the terms at each power of 5. 
As the result one obtains the set of equations 

d 2 x + q (x) - A(fc)] * n (s, k) = R n (x, k), (19) 

where the inhomogeneous term R n {x,k) can be expressed in terms of functions *&i(x, k) 
at i < n — 1, Xi(k) and \ii at % < n and the potential qi(x). For instance, the first two 
terms have the form 

Ri(x, k) = [Xx(k) - qi(x) + 2im(k)d x ] $> (x, k), 
R 2 (x,k) = [Ai(fe) -qi(x) + 2im{k)d x ] ® 1 {x,k)+ 



+ 



X 2 {k)-n\{k) + 2in 2 {k)d x ^ (x,k) 



Equation ((T^)) is to be solved in the class of functions with quasiperiodic boundary condi- 
tions 

V n (x + c,k) = e ikc V n (x,k). (20) 

According to the Fredholm theorem the solvability conditions of eq. (|19|) is the orthog- 
onality of the r.h.s. of (|19|) to the solutions of the corresponding homogeneous equation 

-d 2 x +q (x)-X(k)]¥^(x,k)=0 (21) 

with the same boundary conditions (|2"U|) . It is obvious that the unique solution of eq. (fTU|) 
satisfying the boundary conditions (|2()|) is given by the Bloch function ^o(^) k). Thus the 
solvability conditions have the form 

c 

(R n , * ) = J Rn(x, k)$ (x, k)dx = 0. (22) 



The n-th solvability condition from (1221) gives the second linear link between X n {k) and 
H n {k). 

Before writing down it explicitly let us note that the solvability conditions of equations 
similar to (|19|) has been used in [01 ll()j for constructing the quasiclassical series in the 
problem of the adiabatic perturbation of the periodic potential. In the present paper we 
construct the perturbational series, however the specific character of perturbation leads to 
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the necessity of using the technical tools which are more appropriate for the quasiclassical 
approach. 

Turn back to the link between \ n (k) and fJ> n (k) and introduce the following functions 
of k: 



N n (k)= J ^ n (x,k)^ (x,k)dx, @ n {k) = J qi(x)ty n (x,k)$> (x,k)dx, 



c 

v n {k) = -i J ^' n (x,k)^ (x,k)dx. 


The constructions described above can be summarized in the following 
Proposition 3. The links between X n (k) and fJ, n (k) have the form 

X n (k)N (k)-2v (k)iJ n (k)=f n {X tl ^N t (k),v t (k)(i<n-l),p n ^(k)}. (23) 

where T n are some nonlinear functions of its arguments which can be calculated explicitly 
for arbitrary n. In particular, the first two functions have the form 

fx = Po(k), 

f 2 = 0x(k) + n\N - X X N X + 2ii x vx. 

Remark . To calculate T n for arbitrary n one has to insert R n (x,k) which are known 
explicitly into the solvability condition (|22|). 

Eas. ()18j) and ((2*3*|) give the 2x2 linear system for the functions X n (k) and (J, n {k). By 
solving this system one obtains the following formulae. 

LEMMA 1. The correction terms of expansions hl!3\) and for the quasimomentum 
p and dispersion function E(p) have the form 

= f n X(k) + 2 V0 (k)T 
X(k)No(k)-2v (kY 

u (k) = f " + iVo ^ T " (25) 
X(k)N (k)-2vo(kY 

Here T n and T n are described in propositions 1 and 2 respectively. 

Finally one has to find the correction terms in the series by solving the boundary 
value problem (|19|) . I|20j) . Its solution can be represented as the linear combination of a 
general solution of (|21j) and a particular solution of (|19|) . Thus 

x x 1 

V n (x, k) = B n (k)^ (x, k) - * (x, k) J % 2 dx' I R n ^> Q dx" (26) 

where B n (k), x n and x' n are some constants w.r.t. x. These constants can be defined 
unambiguously by means of quasiperiodic boundary conditions (j2U|) and the normalization 
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condition ((TB|) , It gives x n = 0, B n (k) = and the constant x' n is denned by the equation 

r dx'% 2 r; n R n * dx" = o. 

The expressions (|24|) . (|25j) and (|26|) determine the expansions pip . (|12|) . com- 
pletely. However the correction terms of these expansions have singularities at some values 
of unperturbed quasimomentum k s . As it has been said already we call the corresponding 
points X s (k) as the singular points. 

LEMMA 2. There are two sets of singular points. One of them {Ef}i>i coincides 
with the ends of bands of the unperturbed operator Hq. The singular points {e^i °f 
another set are located inside bands of Hq. If the ratio of the perturbed potential period to 
the unperturbed potential period is m = c/a then inside the N-th band there are (m — 1) 
singular points distributed in such a way that the corresponding values of quasimomentum 
divide the Brilloin zone by m equal segments, so 

$N = *(—), l = l,-,m-l. (27) 
V ma J 

Proof. Firstly let us show that at the ends of bands of the unperturbed operator Hq 
the function Vo(k) vanishes: vo(fcjv) = 0, (fcjv = nN/a). Indeed at the ends of bands the 
Bloch functions *$>Q(x,k) and ^>Q(x,k) coincide: = ^o- Hence 

c c 

v = -ij y'o^odx = -ij y'o^odx = ~ 



As it has been mentioned already at the ends of bands X(k[) = 0. Hence the denominator 
of fractions (|24|) , (|25|) vanishes at the ends of bands and the perturbation series 1)11(1 , (|12|) 
and 1)13(1 are not valid at these points. 

The second set of singular points is generated by the singularities of T n which can 
be shown to be the same as of T2 given by (20) and hence are defined by the equation 
F'(X) = ^F(X) = 0. Since F(X) = coskc one has thatF'(A) = -sin/cc^. and hence 
F'(X) = if sin kc = or dk/dX = 0. In the first case k = ki]y = irlN/ma, I = 
l,...,m, N G Z. However if I = m the quasimomentum k m N corresponds to the end of 
band of the unperturbed operator Hq and the zero of sin k m NC is compensated by the zero 
of A(fc m jv). So we have to consider only kuy with I = l,...,m — 1. These points divide 
each Brilloin zone by m equal segments. The corresponding energy points e® N = X(kipf), 
I = 1, ...,m — 1, are located inside the N-th band. 

The points where dk/dX = lie inside gaps. However the singularities of T n at these 
points are compensated by the factor X(k) entering the fractions and ()25[) . Q.E.D. 

To complete the study of expansions (fi"2|) and ([T3|l outside the singular points 

one needs to analyze the properties of correction terms X n (k) and /J, n (k) when k is inside 
or outside the bands of the unperturbed operator Hq. Remind that to the bands there 
correspond real values of quasimomentum k, while in gaps quasimomentum has nontrivial 
imaginary part: k = 7rN/a + 27 (N is the gap number). 

LEMMA 3. At real k (k / irN/a, N = 1,2,...) the corrections terms X n (k) and 
Hn(k) are real. If k = irN/a + i^ the correction terms X n (k) of the dispersion function are 
real while the correction terms fJ> n (k) of quasimomentum are pure imaginary. 



Jli-nlm ,t,2 



^(0,^)-*o(0,^) 



s 



Remark . The statement of lemma means that outside the singular points the spectral 
bands and gaps of the unperturbed operator Hq are stable under the perturbation. It 
means that if some point A (A / ef, A / ^in) De l° n g s to the band (gap) of Hq the shifted 
point E = A + 0(6) should also belong to the band (gap) of the perturbed operator H. 
Indeed, if k is real then /J, n (k) are also real and the perturbed quasimomentum p given 
by expansion is also real. If k = nN/a + then fi n (k) are pure imaginary and 
p = TrN/a + 1(^ + 0(5)). Hence the corresponding spectral point E(p) belongs to the gap 
of H. It should be noted that the number of this gap is equal to mN since the width of 
the Brilloin zone of the perturbed operator is ir/c = n/ma. 

Proof. Here we prove the statement of lemma for n = 1. For n > 2 the statement can 
be proven analogously. Remind that 



X(k)N (k) -2v (k) X(k)N (k) - 2v (k) 

To prove lemma we show: 1) that at real k the quantities X(k), (3o(k), No(k) and vo(k) 
are real; 2) if k = irN/a + 17 the functions No(k) and /3o(k) are real while the functions 
A(&) and vo(k) are pure imaginary. 

Since cos ka = Fq(X) where Fo(A) is the Lyapounov functions of the unperturbed 
operator Hq, one has 

; d , . , , a sin ka . , , . d „ . , . 

Since the function .F'o(A) is always real at A S R, the function s'mka is real at real k 
and s'm(Nir + ija) = (—1) N isinhja is pure imaginary at k = nN/a + ij, then X(k) is real 
inside bands and is pure imaginary inside gaps. 

The fact that functions No(k) and Po(k) are always real directly follows from their 
definitions. Inside gaps (k = irN/a + £7) the Bloch functions are real. Therefore the 
function 



v (k) = -i I ^'o(x,k)^ (x,k)dx 

is pure imaginary. At real k 







v (k) =i J ^' Q (x,k)^ (x,k)dx = i^' ^ \^-i J <S' 9 dx 




-i J VoVodx = v (k), 




so the function vo(k) is real. Q.E.D. 
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3 Local perturbation series at singular points 



In this Section we construct the local perturbation series at singular points and by the 
study of their properties prove two basic statements. To formulate the first statement we 
need the following notations. Let m = c/a be the ratio of the perturbed potential period 
to the unperturbed potential one. We introduce the points km = irlN/ma, I = 1, m — 1, 
where N is arbitrary integer. According to lemma 2 the points e9 N = X(km) are singular. 
Let us introduce the functions 



and 



<p (x) = - [V (x, km) + # (z, k m )] (28) 

Mx) = [*o(a:, h N ) - V (x, k m )] , (29) 

where ^q{x, k) is the Bloch function of the unperturbed operator Hq fixed by the conditions 
l|15|). Consider the quantities em and pin given up to the second order of 5 as follows 

em = e m + 5ej N + 0{5 2 ) (30) 

and 

Pm = k lN + tift N + 0(S 2 ) (31) 

Here 



1 _ J WWf" fnr,\ 

£ IN — re TYTZXIZ 



i _ H li{xm{x) dx 
Jo <Po( x ) dx 
and 

a\ _ - £ inJo (po(x)(p dx + Jo qi(x)(po(x)<p Q dx 

~ 2/S^(x)dx ■ (33) 

Under these notations the following statement is valid. 

THEOREM l.At the singular points e^, Z = 1, ...,m — 1 f/ie iV-tft spectral band of 
the unperturbed operator Hq is splitted under the perturbation 5qi(x) into m bands. The 
separating gaps are concentrated around the points em given by 13 The quasimomenta 
Pin corresponding to these points are complex and are described by The width of the 

IN-th gap in the leading order of 5 is 

A lN ~ 25\(3} N \ + 0(S 2 ) (34) 

with (3} N given by eq. i.V.V)) . 

Proof. To proof the statement we construct the local perturbation series at the singular 
point e® N . Denote by em the point obtained from e® N under the perturbation 5qi(x). Let 
us seek this point in the form 

£iN = e m + J2 Sn£ iN- (35) 

n>l 
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The corresponding value of quasimomentum p we denote by pim and also seek in the 
form of series 

PIN = kN + Y, (36) 

n>l 

where kiN = -rrlN/ma and e® N = \(kiw). 

According to eq. ((SJ) the function $>(x,p) related with the Bloch function ty(x,p) by 
eq. © at the point p = pin has the property 

*{x + c, PlN ) = e ik ^(x, PlN ) = (-l) lN $(x, PlN ). (37) 

Thus $(x) = $(x,pin) is periodic or antiperiodic function. Let us seek it in the form 

$(x) = <^o(x) + ^5 n ^(x), (38) 

n>l 

where the function ipo (x) is given by (|28jl. Note that due to (|T5|) the leading term (fo(x) 
is periodic or antiperiodic, i.e. (po(x + c) = (— 1) cpo(x). To provide the property (J37f) the 
correction terms ip n {x) are to be sought in a class of functions with the same boundary 
condition: ip n (x + c) = (— l) lN (p n (x). 

After inserting the Bloch function ^(x,pin) = exp{i(piN — kiN)x}Q(x) into equation 
with E 1 = one obtains for the following equation 

-<9^ - 2i(pijv - kiN)9 x + (pin - kiN? + qo(x) + - e [N = 0. 

Now we insert the expansions (|35|) . (|3fij). (|38j) into the last equation and equate terms at 
equal powers of 5. In the leading order of 8 one obtains identity 

-dl + Qo(x) - em ip (x) = 0. 

In the first order one has the equation 

-dl + q (x) + -e m ipi(x) = ej N - q x (x) + 2ijj,j N d x ip (x) = Rx(x). 

Following the discussions of Sec. 2 one concludes that the condition of solvability of this 
equation in a class of periodic (antiperiodic) functions is the orthogonality of R\ (x) to pe- 
riodic (antiperiodic) solutions of the corresponding homogeneous equation. Since there are 
two such solutions (say, for example, ipo(x) and <p(x) given by Q28JI and ((2^|) respectively) 
one has two solvability conditions 



Ri(x)ipo(x)dx = 0, 



Ri(x)ipo(x)dx = 0, 



where the fact that (po and cpo are real is taken into account. Since Jq <p' Q (x)<po(x)dx = 0, 
the first condition allows to define e\ N and yields relation (|32|) . The second condition leads 
to the relation fj,j N = i[5} N with (5} N defined by (|33jl. 
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Now one knows the r.h.s. R\{x) explicitly and can easily find the periodic (antiperiodic) 
solution ip\{x). Then the whole procedure can be repeated. Namely, one can write down 
the equation in the second order of 5 (equation for 992 (x)) and find the correction terms 
ef N and pf N from two solvability conditions 



R2(x)(po(x)dx = 0, J R2(x)Lfo(x)dx = 0. 



The described procedure allows to define the correction terms in expansions , (|3f))) , (jSBj) 
of arbitrary order. However this is out the scope of the present paper. 

Since the width of the Brilloin zone for the perturbed operator is ir/c = n/ma then 
according to the general theory of periodic operators the complex quasimomentum pi^ = 
irlN/ma + iS[3f N + 0(5 2 ) corresponds to the real energy e^y = E(j>in) which lies in the 
IN-th gap of H. 

To calculate the width of gaps appeared under the perturbation let us calculate the 
Lyapounov function T{E) of the perturbed operator at the point E = ejjy: 

Hew) = cos(k lN + i6$ N + 0(S 2 )) = (-l) lN cosh(6\$ N \ + 0{5 2 )) = 
= { _ 1) w ( 1+ S 3M + (S i )). 



Now let us expand the Lyapounov function in the vicinity of euy in the Taylor series 
F(E) = T(e lN ) + T " {£ 2 in) (E - ei N ) 2 + 0((E - e lN ) 3 ) and set E = Ef N where E± N are 
the ends of the gap. Here we used the fact that J-'{£in) = 0. Since the number of the 
considered gap is IN one has 

HEf N ) = F"{e lN ) = (-l) w + 1 |^(e| JV )|. 

and hence the width of the IN-ih gap in the leading order of 5 is 

A m = (E+ N - e w ) - (Ef N - e w ) ~ 25\$ N \. 

Q.E.D. 

Now it remains to consider the singular points which coincide with ends of bands 
(gaps) of the unperturbed operator. Let us denote by E^\. = \{nN/a + 0) and E^_ = 
X(7rN/a — 0) the right and left ends of the iV-th gap in the spectrum of the unperturbed 
operator Hq. The following statement is valid. 

LEMMA 4. Under the perturbation 5q\{x) the ends of the N-th gap of the unperturbed 
operator are shifted as follows 

E in ± = eP ± + 5E%{+0(5 2 ) (39) 



where 

_ S o qi{x)^>l(x )dx 
S c * 2 (x)d. 



h N± - — witnn^iz — \ m ) 
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and *&o(x) = ^o(x, nN/a =L 0) is the Block function ^o(x,k) of Hq taken at the ends of 
the N-th gap. 

Proof. To obtain the expansion (|3*9*)l let us insert the series 

%± = E$l + E = *o(*) + E ^n(a?) (41) 

n>l n>l 

into equation (|1U|) considered at E" = E\n± where Ei?q± are the boundary points of the 
perturbed operator spectrum. Here *$>(x) is the Bloch function of H at the point E[n±. 
In the leading order of 5 one obtains the identity 

-dl + q (x)-E$l] M^)=0. 

The first order term has the form 



-d 2 x + q (x) - E§1\ fa(x) = E^ {x) - q x {x)^{x) = R x . (42) 

Since the functions ^f(x) and ipo(x) are periodic (antiperiodic) with the period c, one 
has to find the correction terms ^ n (x) in the same class. The solvability condition of 
equation (|42j) in such class of functions is the orthogonality of R2 to periodic (antiperiodic) 
solutions of homogeneous equation. Since at the boundary points of spectrum ^o(x, k) = 
^o( x , k), the function 1 i'o(x) is the single periodic (antiperiodic) solution and the solvability 

condition has the form / R2(x)^o(x)dx = 0. It yields the first correction term E N ± ^ n the 



form P0|). 

Obviously the above described procedure can be repeated again and allows to define 
all terms E^}_ in expansion 1)41)1 . Q.E.D. 

It is convenient to combine the statements of theorem 1 and lemma 4 in the following 
form. 

THEOREM 2. Let c/a = m £ N be the ratio of the perturbed operator H period 
to the unperturbed operator Hq one. Then the gaps in the spectrum of H can be divided 
into two sets: the "old" gaps with numbers mN, N = 1,2, ... which are just shifted gaps 
of the unperturbed operator Hq with the shifts given by eas. \3y\) , \4(fij ; and "new" gaps with 
numbers IN, I = 1,2,..., m — 1,(N = 1,2,...) which arise under perturbation. These new 
gaps are located inside the N-th band of the unperturbed operator, are concentrated around 
the points ein given by eg. h30\) and have the widths described by eg. \34\l- 
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